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We present progress in the QCD analysis of the Bjorken sum rule at low momentum transfers. We 
study asymptotic structure of the perturbative QCD expansion at low Q 2 scales based on analysis 
of recent accurate data on the Bjorken sum rule and available now four-loop expression for the 
coefficient-function Cbj(Q 2 )- We demonstrate that the standard perturbative series for Csj(Q 2 ) 
gives a hint to its asymptotic nature manifesting itself in the region Q 2 < 1 GeV 2 . It is confirmed 
by the considered integral model for the perturbative QCD correction. We extract a value of higher- 
twist p,4 coefficient and study the interplay between higher orders and higher-twist contributions. 

£f~^ , Results of other approaches to the description of Bjorken sum rule data are discussed. 
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I. INTRODUCTION 



Perturbation theory (PT), supplemented by the renormalization group method and the renormalization procedure, 
is a main tool of QCD calculations. Conventional perturbative expansions for hadronic observables are the power 
r- — i series in the strong coupling a s . Owing to the property of asymptotic freedom in QCD, the perturbative description 
'■^ ■ becomes more reliable at high energy region. On the contrary, at low momentum transfer, Q 2 < 1 — 2 GeV 2 , the 
use of the perturbative power series may be questionable. Therefore, it is important to know at least a few terms of 
I- PT series to estimate theoretical uncertainties associated with a truncation of the perturbative expansion. Currently, 
the perturbative QCD analysis is aimed at the perturbative calculations of new terms of PT series. Only recently up 
to the fourth-order perturbative approximation for some physical processes became available 0, H| . The number of 
these processes includes the fundamental sum rule for polarized deep inelastic scattering which is the Bjorken sum 
t-H \ rule Pi _ "(Q 2 ) [3j defined as a difference between the first moments of the proton g\ and the neutron spin structure 
functions. The Bjorken sum rule has been measured in polarized deep inelastic lepton scattering at SLAC, CERN, 
DESY At low Q 2 , in a Q 2 -range from 0.05 < Q 2 < 3.2 GeV 2 , high-precision data on the T\~ n has been 

presented by the Jefferson Lab [l0[ and at 3.0 GeV 2 by the COMPASS collaboration [ll|. 



in 

At low momentum scales, the QCD analysis of the Bjorken sum rule includes both the perturbative and non- 
perturbative higher-twist (HT) components related to each other. It is clear that the reliability of extracting infor- 
mation on the HT effects is connected to the indeterminacy in the description of the PT series. As new perturbative 
correction of an order of to the Bjorken sum rule became available 0, the opportunity for new researches opened. 

In this report, we present results of the four-loo p Q CD analysis of the Bjorken sum rule at low momentum scales, 
continuing our investigations started in Rcfs. [12, fl3j | - We consider the features of the four- loop PT series and the 
interplay between the higher PT order and the HT contributions. We discuss results of application of other approaches 
which were applied to the description of low Q 2 Jefferson Lab data. 



II. STANDARD APPROACH 



Away from the large Q 2 limit, the polarized Bjorken sum rule is given by a double series in powers of a s and in 
powers of l/Q 2 (nonperturbative HT corrections): 



r p - n (Q 2 ) = M 



i=2 



OQ 

^ M2i (I) 



where \g^\ = 1.2701 ± 0.0025 14[ is the nucleon axial charge, /X4, . . . are the HT coefficients; A^J(Q 2 ) is the 
perturbative correction which is denned by the coefficient-function CsjiQ 2 ), A^J(Q 2 ) = 1 — Cej (Q 2 )- Note that until 
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very recently A^J(Q 2 ) has been known up to the third order in a s . The corresponding expression was used in many 
studies (see, e.g., Refs. [H, hi particular, to extract a value of a s from experimental data. Comparison of 

these values with other accurate a s values, such as those obtained from the r-lepton and the Z-boson into hadrons 
width decays, is an important test of the consistency of QCD [lil. Il8j. 



A. Four-loop analysis 
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FIG. 1: The relative contribution N;(<2 2 ) = 
Cial(Q 2 )/ABj(Q' 2 ) for the i-th term of series © as 
a function of the Q 2 . 



FIG. 2: Comparison of the 1-, 2-, 3-, and 4-term approxi- 
mants of series ((4]) with the exact result <(3j and with N 3 LO 
Bjorken sum rule series ([2]). 



At the four-loop (N 3 LO) level, the perturbative QCD correction to the Bjorken sum rule in the case of massless 
quarks in the MS rcnormalization scheme for / = 3 light quarks flavors reads as Q 

A^ T = 0.318 a s + 0.363 a 2 + 0.652 a 3 + 1.804 a* . (2) 

We use the PT running coup ling a s (Q 2 ) obtained by integration of the renormalization group equation with the 
four- loop (3- function (see Ref. [13| for additional details). 

Let us discuss the convergence properties of the PT power series ^ at low momentum transfers. Figure [T] shows 
the relative contribution of the i-th term of this series as a function of Q 2 . As can be seen from this figure, in the 
region of small Q 2 < 1 GeV 2 the dominant contribution comes from the a 4 -term. This may be considered as an 
indication of the transition of PT series to the asymptotic regime and one can estimate the value a s (l GeV 2 ) ~ 0.5 as 
a critical. In the region Q 2 > 2 GeV 2 the situation changes - the major contribution comes from the first and second 
terms. 

To specify discussed above convergence properties of series ([2]) , we consider the exact example following from the 
integral model [l9| (see Ref. [2(| for more details): 



4 fe r(fc + 1) 



Using this expression, one can build series which is close to the series 0: 



bu 

A(g) = ^- [C(g) - 1] = 0.318 g + 0.348 g 2 + 0.718 g 3 + 2.188 g 4 + . . . . (4) 

37T 

In Fig. [2] we compare the finite sum approximations A[„](<7) = ^2 n Cig l , n — 1, ... ,4 both with the exact result 
and with the Bjorken series ([2]). This figure demonstrates the closeness of the 4-term approximation (dashed curve) 
and the N 3 LO approximation ([2]) (dotted curve). It is necessary to note, that 1-, 2-, and 3-order approximations to 
A(g) and Aej practically coincide with each other. This figure shows that the 2-term approximant is good up only 
to g = 0.15 — 0.20 and the 3-term one up to g ~ 0.33 while the 4-term sum starts to deviate from exact C(<?)-curve 
at g ~ 0.27 . We interpret this observation as asymptotic structure manifestation. This model confirms that the 
asymptotic structure of the pQCD expansions manifests. Therefore, the application of standard PT greater than 
N~LO approximation can not allows to extract accurate information in the low-energy domain. 
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B. Higher twist contribution 



We expand our consideration, including the HT part which is presented in expression ([I} . In Table Q] we show our 
results for values of the coefficient /i 4 (the errors are statistical only) fitted to the low Q 2 data [Hj, EH in different 
PT orders. One can see that /Z4-values extracted changes rather strongly between different PT orders. The absolute 
value of fi,4 decreases with the order of PT and just at N 3 LO becomes compatible to zero. It can be interpreted 
as a manifestation of duality between higher orders and HT (see Ref. Note that a value extracted in the 



TABLE I: Results of /M-extraction from the data on the Bjorken sum rule in different PT orders. 



PT order 


LO 


NLO 


N 2 LO 


N 3 LO 


tn, GeV 2 


- 0.037 ± 0.003 


- 0.025 ±0.004 


-0.012 ±0.006 


0.005 ± 0.008 



leading-order (LO) is consistent with a value /i4 = — 0.047 ± 0.025 GeV 2 presented in Ref. [22J as well as with a value 
/14 = —0.028 ± 0.019 GeV 2 obtained from the next-to-leading-order (NLO) fit based on the ^-dependent structure 
functions data [23[ . 



III. OTHER APPROACHES 



As shown above, use of the conventional PT series even with the HT component does not allow to describe the 
Jefferson Lab data down to the infrared region (see Fig. [3]). Let us consider some other approaches which are more 
successful in this direction. 



A. APT 



In Ref. 2J| the conventional method of the renormalization group improvement of the perturbative expansions was 
modified by requiring Kallen-Lehmann analyticity, which reflects the principle of causality. In the framework of this 
approach !24j called as the analytic perturbation theory (APT) the ghost pole and corresponding branch points, which 
appear in higher PT orders, are absent (see, e.g., Ref. [Hj]). As the moments of the structure functions should be 
analytic functions in the complex Q 2 plane with a cut along the negative real axis (see Ref. [26[ for more details), the 
standard description of the r^ - "(Q 2 ) violates analytic properties due to the unphysical singularities of perturbative 
running coupling. On the other hand, the APT support these analytic properties (see Ref. [27[ as review). The 
perturbative correction to the Bjorken sum rule can be written in the form of a spectral representation [l6j and at 
four-loop level is 

A^f T = 0.318 Ai + 0.363 A 2 ±0.652 _4 3 ± 1-804 At. (5) 

The power coefficients in this expression are the same as in series @, and the functions Ak{Q 2 ) are defined through 
the spectral density Qk(cr) = Im [a* (— a — ie)] by the spectral integral (see Ref. [l3j for additional details). At large 
momentum transfers, all the functions Ak(Q 2 ) become proportional to the fc-th power of the usual perturbative 
coupling [a s (Q 2 )] fc and the expansion reduces to the power series ©. However, at small enough Q < 1 — 2 GeV 
the properties of the non-power expansion ([5]) become considerably different from the PT power series ([2]). 

Figure [3] shows the results of /i4-fits in different PT orders both in the PT and APT approaches. From this figure 
follows that in the framework APT including /14-coefficient allows one to move further down to Q 2 ~ 0.1 GeV 2 [l3[ 
in description of the Jefferson Lab data. It is important to note that the APT leads to higher-loops stability of the 
HT-extraction: [i^ — —0.044 ± 0.004 GeV 2 in all loop approximations. 



B. MPT 

In context of data on the r^ _ ™ analysis it should be noted the MPT approach 28]. Basis of the MPT is a simple 
idea to change the usual logarithm in the expression for the QCD running coupling, 1/ ln(Q 2 /A 2 ) , that is singular in 
the infrared region, on the "long logarithm" ln(£ ± Q 2 /A 2 ) , where parameter £ corresponds to the "effective gluonic 
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mass" m g i = y/£ A: 

A hM PT(Q 2 ) = J o Hc l Q2/A2 y (6) 

Result of MPT application can be interpreted as a transition to the new momentum-transfer scale both in the PT 
and the HT contributions: A^j PT ~ A^ T (Q 2 + m 2 gl ) and n±/Q 2 ~ ^/ (Q 2 + m l t ) ■ Tne MPT curve is shown in 
Fig. 2] as solid line (preliminary result Ref. [II])- From this figure one can see that the MPT approach allows us to 
describe the Bjorken sum rule data down the infrared region Q 2 — > 0. 

Note that very close description of the Jefferson Lab data was obtained in Ref. [23], where was used a modified 
expression for the coefficient- function Cej in a combination with "frozen" behavior of running coupling a s {Q 2 ). 



IV. SUMMARY 



We have presented the QCD analysis of the Bjorken sum rule up to four-loop level in a s at low momentum transfers. 
It has been shown that the asymptotic nature of the perturbative series manifested itself at the four-loop level in 
the region Q 2 < 1 GcV 2 . Therefore, at low Q 2 the inclusion of the four-loop term does not improve the precision of 
theoretical predictions. It was confirmed by the considered integral model for the perturbative QCD correction. 

Using the recent low Q 2 -data from the Jefferson Lab and COMPASS experiments, we have extracted a value of 
higher-twist \x± coefficient and have shown the interplay between higher orders and higher-twist contributions. Results 
of other approaches to the description of Bjorken sum rule data have discussed. 
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